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A  Plume  Comparison  of  Xenon  and  Krypton 
Propellant  on  a  600  W  Hall  Thruster 

Michael  R.  Nakles*  Ryne  R.  Barry! 

ERC,  Inc.,  Edwards  Air  Force  Base,  CA  93524 

C.  William  Larson*  and  William  A.  Hargus,  Jr.* 

Air  Force  Research  Laboratory,  Edwards  Air  Force  Base,  CA  93524 


An  experimental  study  was  conducted  to  compare  variations  in  the  plume  of  a  600  W 
Hall  effect  thruster  operating  on  xenon  and  krypton  propellants.  Both  a  guarded  Faraday 
probe  and  a  retarding  potential  analyzer  were  used  to  measure  plume  properties  for  vari¬ 
ous  thruster  operating  conditions.  Discharge  current  data  was  also  measured  in  the  time 
domain  using  an  oscilloscope  and  in  the  frequency  domain  using  a  signal  analyzer  for  all 
tested  operating  conditions.  Operating  conditions  were  chosen  to  provide  a  comparison  of 
the  effects  of  radial  magnetic  field  strength  for  both  propellants.  For  each  xenon  propel¬ 
lant  operating  condition,  two  krypton  propellant  operating  conditions  were  tested.  One 
condition  matched  volumetric  flow  rate  and  the  other  condition  matched  anode  power  by 
increasing  the  anode  flow  rate.  The  plume  data  were  analyzed  to  calculate  various  com¬ 
ponents  of  the  energy  and  propellant  utilization  efficiencies  as  defined  in  Ref.  1.  Faraday 
probe  data  were  used  to  calculate  the  divergence  and  beam  current  utilization  efficiencies. 
Voltage  utilization  efficiency  was  calculated  using  data  from  the  retarding  potential  ana¬ 
lyzer.  Compared  with  xenon,  krypton  divergence  efficiency  is  approximately  15-20%  lower 
and  its  current  utilization  efficiency  is  5-10%  lower.  Differences  between  voltage  utilization 
efficiency  for  xenon  and  krypton  were  small.  Anode  current  oscillation  magnitude  was 
lower  for  krypton  propellant. 


Nomenclature 


e  elementary  charge 

Ia  anode  discharge  current 

A  integrated  thruster  beam  current 

Im i  magnet  current  to  the  four  series  outer  magnetic  cores 
Im2  magnet  current  to  the  central  magnetic  core 
j  plume  charge  flux 

m  atomic  mass 

ma  anode  mass  flow  rate 

P  anode  power 

q  multiple  for  units  of  elementary  charge 

r  radial  distance  in  thruster  coordinate  system 

V  ion  acceleration  voltage 

V  mean  ion  acceleration  voltage  from  a  distribution  function 

Va  anode  potential  measured  with  respect  to  cathode  common  potential 
Va  anode  volumetric  flow  rate 

Vo  ion  retarding  grid  potential  for  retarding  potential  analyzer 

‘Research  Engineer,  AFRL/RZSS,  1  Ara  Rd.  Edwards  AFB,  CA  93524 
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r  current  collected  by  retarding  potential  analyzer 
AV  average  ion  acceleration  voltage 
r]a  anode  efficiency 

6  angular  position  in  thruster  coordinate  system 

<  e  >j  charge  flux  weighted  plume  divergence  angle 

*dIV  divergence  loss  component  of  propellant  utilization 

<i>vdf  velocity  distribution  loss  component  of  propellant  utilization 

<>j  charge  flux  weighted  average  quantity  in  the  plume 

<>m  mass  weighted  average  quantity  in  the  plume 

<>mv  momentum  weighted  average  quantity  in  the  plume 


Introduction 

Due  to  a  number  of  engineering  reasons,  xenon  (Xe)  is  the  propellant  of  choice  for  Hall  effect  thrusters. 

These  include  its  high  mass  (131  amu)  and  its  relatively  low  ionization  potential  (12.1  eV).  Furthermore, 
the  inert  nature  of  xenon  eliminates  much  of  the  controversy  that  plagued  early  electrostatic  propulsion 
efforts  when  mercury  (Hg)  and  cesium  (Cs)  were  the  propellants  of  choice.  Although  xenon  is  a  noble 
gas,  it  is  the  most  massive,  and  due  to  its  non-ideal  gas  behavior,  it  is  possible  to  pressurize  and  store  at 
specific  densities  that  exceeded  unity.  As  such,  it  can  be  stored  at  higher  densities  than  the  common  liquid 
monopropellant  hydrazine.2 

While  xenon  will  likely  remain  the  ideal  propellant  for  electrostatic  electric  propulsion  thrusters.  There 
are  several  concerns  that  have  driven  the  Hall  effect  thruster  community  to  explore  alternative  propellants. 
As  orbit  raising  missions  of  longer  duration  and  larger  payloads  are  proposed  for  Hall  effect  thrusters,  the 
mass  of  required  propellant  increases.  Xenon  production  is  a  byproduct  of  the  fractional  distillation  of 
atmospheric  gases  for  use  primarily  by  the  steel  industry.  Due  to  the  low  concentration  of  xenon  in  the 
atmosphere  (~  90  ppb),  worldwide  production  appears  to  be  limited  to  approximately  6,000  m3  per  year. 
Increasing  industrial  demand  for  items  such  as  high  efficiency  lighting  and  windows  has  produced  wide  price 
swings  in  the  past  decade.  Xenon  prices  have  varied  by  as  much  as  a  factor  of  ten. 

For  high  thrust  to  power  missions,  bismuth  (Bi)  has  been  demonstrated  as  a  viable  alternative  Hall  effect 
thruster  propellant.  Bismuth,  with  its  high  atomic  mass  (209  amu)  and  low  ionization  potential  (7.3  eV) 
appears  to  have  advantages  for  missions  where  high  thrust  at  reduced  specific  impulse  is  advantageous, 
primarily  orbit  raising  missions.  Bismuth’s  main  drawback  is  that  the  metal  must  be  vaporized  to  be 
ionized  and  accelerated  within  a  Hall  effect  thruster.  The  requirement  for  high  temperatures  (boiling  point 
of  1,837  K)  require  special  engineering  considerations  compared  to  the  relatively  simple  gas  distribution 
systems  used  for  xenon.  In  addition,  the  use  of  vapor  as  a  propellant  has  tended  to  cause  concern  for 
spacecraft  operators  despite  the  assurances  of  thruster  developers. 

For  missions  that  can  benefit  from  higher  specific  impulse,  krypton  (Kr)  may  have  some  benefits.  Krypton 
has  a  lower  atomic  mass  (83.8  amu),  but  a  higher  ionization  potential  (14.0  eV)  than  xenon.  However  like 
xenon,  krypton  is  a  noble  gas  and  could  be  easily  integrated  into  existing  xenon  propellant  management 
systems  without  much  modification.  The  similar  ionization  potential  would  likely  not  dramatically  affect 
the  efficiency  of  a  Hall  effect  thruster,  and  the  lower  mass  would  produce  a  25%  increase  in  specific  impulse 
assuming  there  were  no  offsetting  losses.  The  increase  in  specific  impulse  would  be  useful  for  missions  such 
as  station-keeping.  For  missions  such  as  orbit  raising,  increasing  the  specific  impulse  will  increase  trip  time 
due  to  power  limitations.  However  as  solar  electric  power  system  specific  power  decreases,  increasing  the 
specific  impulse  of  the  propulsion  system  is  advantageous.  Krypton  is  approximately  10  times  more  common 
in  the  atmosphere  (and  hence  in  production)  than  xenon,  and  when  accounting  for  mass  is  approximately  6 
times  less  expensive.  One  disadvantage  for  krypton  is  that  its  tankage  fraction  appears  to  be  substantially 
higher  than  that  of  xenon  due  to  reduced  van  der  Waals  interactions.  As  such,  compressed  gas  tankage 
fractions  could  be  as  high  as  37%.  At  least  one  study  has  examined  this  issue  and  has  identified  space  rated 
cryo-coolers  that  could  liquefy  krypton  (120  K  boiling  point),  or  for  that  matter  xenon  (165  K  boiling  point), 
and  reduce  tankage  fractions  to  less  than  2%. 2 

The  goal  of  this  study  is  to  understand  the  differences  between  xenon  and  krypton  plumes  in  a  medium 
power  Hall  effect  thruster  in  a  design  optimized  for  xenon.  Several  plume  properties  are  measured  and 
some  generalizations  of  various  performance  parameters  are  inferred.  However,  this  study  is  primarily  an 
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exploratory  investigation  into  the  possibilities  of  using  krypton  as  a  propellant  in  Hall  effect  thrusters, 
primarily  due  to  its  low  cost  and  possible  application  as  a  low-development  overhead  replacement  for  existing 
xenon  Hall  effect  thruster  systems. 

Experimental  Apparatus  and  Techniques 


Test  Facility 

The  tests  performed  in  this  study  utilized  Chamber  1  at  the  Air  Force  Research  Laboratory  at  Edwards 
AFB.  Chamber  1  is  a  cylindrical  non-magnetic  stainless  steel  vacuum  chamber  with  a  2.4  m  diameter  and 
a  4.1  m  length.  Pumping  is  provided  by  two  liquid  nitrogen  baffled  (70  K),  1.2  m  flanged  gaseous  helium 
two  stage  cryogenic  (15  K)  vacuum  pumps  with  a  measured  pumping  speed  on  both  xenon  and  krypton  of 
48,500  L/s.  Chamber  pressure  is  monitored  with  a  cold  cathode  gauge.  Background  pressure  for  thruster 
operation  was  measured  to  be  9  x _  6  Torr  for  xenon  and  7  x“  6  Torr  for  krypton  with  gas  correction  factors 
applied. 

The  interior  of  the  chamber  is  covered  with  nuclear  grade,  low  sulfur,  flexible  graphite  4.7  mm  thick.  Both 
chamber  ends  contain  louvered  beam  dumps  manufactured  in-house  using  13  mm  thick,  15  cm  wide  graphite 
panels  to  reduce  redeposition  of  sputtered  materials  on  the  thruster  during  extended  firings.  The  chamber 
floor  is  protected  using  a  carbon-carbon  woven  blanket  that  allows  for  ease  of  placement  and  removal. 

The  thruster  is  mounted  on  a  1.2  m  long  by  1.8  m  wide  test  cart  that  is  capable  of  traversing  on  internal 
chamber  rails  or  being  removed  from  the  chamber  using  a  roller  system.  For  these  tests,  the  thruster  is 
mounted  atop  a  column  and  cantilevered  30  cm  along  the  chamber  central  axis.  A  stepper  motor  driven 
rotary  stage  is  placed  with  its  center  of  rotation  directly  below  the  thruster  exit  plane  directly  on  the  working 
surface  of  the  test  cart.  Due  to  the  column  and  cantilever  placement  of  the  thruster,  the  arm  mounted  on 
the  rotary  stage  is  capable  of  rotating  between  ±150°  from  the  plume  axis. 

The  arm  mounted  on  the  rotary  stage  is  a  stepper  motor  driven  linear  translational  stage  capable  of 
traversing  up  to  80  cm  from  the  center  of  rotation  of  the  rotary  stage.  This  combination  of  rotary  and  linear 
motion  allows  for  the  continuous  sweeping  of  various  plasma  interrogation  probes  in  the  r  —  8  plane. 


(a)  Photograph  of  the  testing  environment  inside  of  (b)  The  probe  translation  system  is  composed  of  a 
Chamber  1.  rotary  stage  for  angular  motion  and  a  linear  stage  for 

radial  movement. 

Figure  1.  Probe  measurements  were  conducted  in  Chamber  1  using  a  stepper  motor-driven  translation  stage 
system  for  movement. 

A  computer  data  acquisition  system  monitors  voltages  and  currents  of  the  power  supplies  used  in  Hall 
thruster  operation.  Data  was  sampled  at  a  rate  of  2  Hz  for  this  study. 

Hall  Effect  Thruster 

The  Hall  thruster  used  in  this  study  was  the  rectangular  600  W  Busek  Company  BHT-HD-600  Hall  thruster 
with  a  3.2  mm  hollow  cathode.  Photographs  of  the  thruster  are  shown  in  Fig.  2.  This  thruster  has  a  conven¬ 
tional  five  magnetic  core  (one  inner,  four  outer)  magnetic  circuit.  The  acceleration  channel  of  the  thruster 
has  a  32  mm  outer  radius  and  a  24  mm  inner  radius  and  a  depth  of  approximately  10  mm  between  the 

3  of  21 


The  31th  International  Electric  Propulsion  Conference,  Ann  Arbor,  Michigan,  USA,  September  20-24,  2009 
Distribution  A:  Approved  for  public  release;  distribution  unlimited. 


Preprint 


geometrical  exit  plane  and  the  furthest  forward  extent  of  the  anode.  The  thruster  has  been  extensively  char¬ 
acterized  to  have  a  thrust  of  39  mN  with  a  specific  impulse  of  1,530  s,  yielding  an  efficiency  of  approximately 
50%  at  the  nominal  conditions  specified  in  Table  1. 


(a)  Busek  BHT-HD-600  Hall  effect  (b)  Thruster  operating  on  krypton  propel- 
thruster.  lant. 


Figure  2.  Photographs  of  the  Busek  BHT-HD-600  Hall  effect  thruster  used  in  this  study. 


Guarded  Faraday  Probe 

Ion  current  flux  was  measured  using  a  guarded  Faraday  probe  pictured  in  Fig.  3.  The  electrodes  were 
constructed  from  molybdenum.  Ion  current  was  collected  with  a  disk  measuring  8.3  mm  in  diameter.  A 
concentric  guard  piece,  measuring  22.5  mm  in  outer  diameter,  was  used  to  minimize  the  effects  of  the  plasma 
sheath  on  the  ion  current  collector’s  effective  collecting  area.  A  0.56  mm  wide  gap  existed  between  the  outer 
wall  of  the  collector  and  the  inner  wall  of  the  guard  ring.  The  effective  current  collector  area  of  the  probe 
was  calculated  by  adding  a  portion  of  the  gap  surface  area  to  the  collector  surface  area.  The  amount  of 
gap  area  that  contributed  to  the  effective  current  collector  area  was  proportional  to  the  ratio  of  lateral  wall 
surface  area  of  the  collector  to  the  total  lateral  wall  surface  area  on  both  sides  of  the  gap  as  suggested  in 
Ref.  3.  Ion  charge  flux  was  measured  by  dividing  the  current  to  the  collector  by  its  effective  surface  area. 
The  disk  and  guard  ring  were  biased  to  -30  V  with  respect  to  chamber  ground  during  the  measurements  so 
that  ion  saturation  was  achieved.  The  effects  of  secondary  electron  emission  were  assumed  to  be  less  than  a 
few  percent4  and  were  neglected  in  the  analysis  of  the  measurements. 


Table  1.  Nominal  BHT-HD-600  Hall  thruster  operating  conditions  for  xenon  propellant. 


Parameter 

Value 

Anode  Flow 

2.45  mg/s  Xe 

Cathode  Flow 

197  f-ig/s  Xe 

Anode  Potential 

300  V 

Anode  Current 

2.16  A 

Keeper  Current 

0.5  A 

Heater  Current 

3.0  A 
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Figure  3.  The  guarded  Faraday  probe  used  in  charge  flux  measurements. 


Retarding  Potential  Analyzer 

The  retarding  potential  analyzer  (RPA)  used  in  this  experimental  investigation  was  a  four  grid  design 
described  in  detail  in  Ref.  3.  The  RPA  operates  as  an  energy  filtered  Faraday  probe.  The  energy  filtering  is 
accomplished  using  a  series  of  biased  grids  in  front  of  the  current  collector. 

The  first  grid  acts  as  an  aperture  to  reduce  plasma  flow  into  probe  and  is  allowed  to  float  at  the  local 
plasma  floating  potential.  The  second  grid  (electron  repulsion),  biased  negatively  with  respect  to  chamber 
ground,  repels  incoming  electrons  within  the  plasma  while  allowing  ions  to  pass  through.  The  potential  of 
the  ion  retarding  grid  (third  grid)  is  swept  positively  relative  to  chamber  ground  to  allow  selective  passage  of 
ions  based  on  their  kinetic  energy.  The  fourth  grid  (electron  suppression)  is  placed  in  front  of  the  collector 
to  suppress  ion  impact  induced  secondary  electron  emission. 

Current  to  the  collector  is  measured  as  a  function  of  ion  retarding  grid  potential.  The  probe  current  is 
differentiated  with  respect  to  ion  retarding  potential  to  calculate  energy-per-charge  distribution  function. 
Analytically,  for  a  single  species,  the  derivative  of  the  current  per  unit  area  is5 


dT 

dVo 


qe 

m 


f(V) 


which  is  proportional  to  the  ion  energy-per-charge  distribution  of  the  plasma. 


(1) 


Electron  Suppression  Grid 
Electron  Repulsion  Grid 
Floating  Grid 


Ion  Retarding  Grid 


(a)  Photograph  of  RPA. 


(b)  Exploded  view  of  the  RPA  showing  grid  placement. 


Figure  4.  Retarding  potential  analyzer  (RPA)  used  to  measure  energy  distribution  functions  (EDF). 


Performance  Evaluation  Techniques 

Reference  1  outlines  a  standard  methodology  for  Hall  thruster  anode  efficiency  analysis.  This  method 
separates  anode  efficiency  into  the  product  of  energy  efficiency  and  propellant  utilization  efficiency.  The 
anode  efficiency  can  be  expressed  as 
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V,  = 


J'VDFY'DIV 


(2) 


where  the  first  two  terms  comprise  the  energy  efficiency  and  the  last  two  terms  produce  the  propellant 
utilization  efficiency.  Energy  efficiency  accounts  for  losses  due  to  Joule  heating,  radiation,  and  ionization 
processes.  It  is  completely  separated  from  jet  vector  properties.  Propellant  utilization  efficiency  accounts 
for  losses  due  to  geometric  plume  divergence  and  ion  species  composition. 

Faraday  probe  measurements  can  be  used  to  approximate  both  the  charge  utilization  efficiency  and  the 
plume  divergence  efficiency.  The  beam  current  can  be  calculated  as 


/h  «  7rr2  j  j(0)|sin#|d0  (3) 

The  ratio  of  beam  current  to  the  anode  current  (recorded  by  the  data  acquisition  system)  is  the  charge 
utilization  efficiency. 

Plume  divergence  efficiency  is  a  measure  of  jet  momentum  loss  due  to  beam  divergence.  It  is  defined  as 
the  ratio  of  axial  momentum  that  produces  thrust  to  the  total  momentum  exhausted  in  the  plume.  It  can 
be  approximated  using  charge  flux-weighted  average  divergence 


=<  cos  8  >1 


7ir2  Jjv  j(8)  cos  8  | sin  8\dd 
nr2  /_%  j{8)  |sin0|d0 


(4) 


Voltage  utilization  efficiency  is  the  ratio  of  average  ion  acceleration  voltage  to  anode  voltage.  Being  a  part 
of  the  energy  efficiency,  it  contains  no  information  about  the  vector  properties  of  the  jet.  Voltage  utilization 
can  be  calculated  as 


/  AV\  _  /<  V{8)  _  l_  / nr2fj„  j(fl)V(fl)|sinfl|dfl\ 

\VaJ  v  Va  )  V„  1  7rr2  J_2k  j{8)  sin  8\d8  J 

where  a  charge  flux-weighted  approximation  enables  calculation  through  a  combination  of  RPA  and  Faraday 
probe  data. 


Results  and  Discussion 


Magnetic  Field  Effects 

Probe  measurements  were  taken  for  four  different  magnetic  field  strength  configurations  for  each  propellant. 
Varying  the  magnetic  field  by  changing  the  applied  current  to  the  four  series  connected  outer  magnetic  cores 
(ml)  and  the  central  core  (m2)  changed  both  the  plume  focusing  and  anode  current  oscillatory  conditions. 
Three  of  these  conditions  have  been  tested  in  the  past  with  xenon. 6-8  The  Im ,  =  4.00  A,  Im2  =  4.00  A 
case  has  the  strongest  radial  magnetic  field  strength.  The  Iml  =  4.00  A,  Im2  =  2.00  A  condition  has  85% 
of  its  magnetic  field  strength  and  the  Im\  =  1.75  A,  Im2  =  1.75  A  case  has  50%  of  the  field  strength. 
Despite  the  change  in  field  strength,  the  shape  of  the  radial  field  profile  remains  similar.  Figure  5  shows 
the  normalized  radial  field  profiles  along  the  acceleration  channel  centerline  calculated  using  a  commercial 
software  simulation  package. 

In  addition  to  the  previously  used  configurations,  an  extra  condition  was  tested  which  attempted  to 
minimize  the  divergence  of  the  beam.  The  magnet  currents  were  tuned  to  maximize  the  current  collected 
by  the  Faraday  probe  at  a  location  of  r  =  80  cm  and  0  =  5°  for  each  propellant.  The  tuned  magnet 
configurations  were  Im i  =  1.07  A,  Im2  =  2.45  A  for  xenon  and  Im ,  =  0.56  A,  Im2  =  2.16  A  for  krypton 
(tuned  in  flow-matched  setting).  (These  configurations  are  referred  to  as  4>div  optimized  in  the  text  and 
figures.) 

For  each  xenon  propellant  operating  condition,  two  corresponding  krypton  propellant  operating  condi¬ 
tions  were  tested  as  in  Ref.  9.  One  condition  matched  anode  volumetric  flow  rate  and  therefore  particle  flux 
through  the  anode.  The  other  condition  matched  anode  power.  Here  the  anode  flow  rate  was  increased  to 
compensate  for  the  smaller  anode  current  per  volumetric  flow  rate  observed  when  operating  with  krypton 
propellant. 
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Figure  5.  Normalized  radial  magnetic  field  strength  along  the  acceleration  channel  centerline.  The  radial  field 
strength  is  a  factor  of  two  greater  in  the  high  strength  case  than  in  the  low  case.  However,  the  profile  of  the 
radial  field  strength  is  nearly  identical  between  both  cases. 


Guarded  Faraday  Probe 

For  most  operating  conditions,  guarded  Faraday  probe  sweeps  were  taken  in  radial  increments  of  10  cm  from 
r  =  20  cm  to  r  =  80  cm.  The  angular  range  of  the  sweeps  was  from  6  =  -120°  to  9  =  120°.  Measurements 
were  taken  with  a  1°  angular  resolution. 

At  each  radius  a  value  for  beam  current  and  beam  divergence  efficiency  were  calculated  using  equations 
3  and  4  integrated  from  -90°  to  90°.  These  calculated  values  were  observed  to  have  a  dependency  on  radius. 
This  dependency  was  stronger  for  beam  divergence  than  total  beam  current.  The  calculated  beam  divergence 
efficiency  value  was  observed  to  decrease  linearly  with  increasing  radius.  This  phenomenon  was  likely  due 
to  collisions  in  the  plume  which  led  to  defocusing.  Integrated  beam  current  slightly  increased  linearly  with 
increasing  radius.  The  greater  population  of  slowly  moving  charge-exchange  ions  at  higher  radii  may  have 
artificially  increased  current  to  the  negatively  biased  probe.  Several  calculated  beam  current  values  for 
xenon  significantly  deviated  from  the  trend  of  the  other  data  points.  Linear  extrapolation  of  the  data  to  the 
cathode  plane  (r  =  2  cm)  was  used  to  arrive  at  representative  values  of  beam  current  and  beam  divergence 
efficiency  for  each  operating  condition.  This  study  differs  from  most  other  studies  in  that  Faraday  probe 
data  was  taken  at  multiple  radii. 

Table  2  summarizes  test  operating  conditions  and  tabulates  the  calculated  values  of  beam  current  and 
divergence  efficiency.  Figures  6,  7,  and  8  show  the  charge  flux  data  and  the  calculated  values  of  beam 
current  and  divergence  efficiency  for  xenon  and  krypton  power-matched  cases.  In  general,  the  integrated 
beam  current  values  were  higher  than  expected.  This  discrepancy  is  believed  to  be  related  to  uncertainty 
in  the  effective  area  of  the  Faraday  probe  current  collector.  This  type  of  systematic  error  would  result  in 
inaccuracy  by  a  consistent  factor  throughout  the  data.  However,  due  to  the  mathematics,  it  would  not 
affect  the  calculation  of  divergence  angle  nor  the  relative  comparison  of  beam  current  and  current  utilization 
efficiency  between  cases. 

Figure  6  shows  Faraday  probe  results  for  xenon.  The  4>div  optimized  magnet  current  setting  and  the 
I  mi  =  1-75  A,  Jm2  =  1.75  A  case  both  resulted  in  the  best  beam  focusing  while  the  Im\  =  4.00  A,  Jm2  = 
2.00  A  case  had  the  most  divergence.  However,  in  beam  current,  the  results  show  a  reversal  in  ranking  among 
the  cases  where  the  4>div  optimized  case  produced  the  lowest  current.  In  terms  of  both  current  utilization 
and  beam  divergence,  the  Im i  =  4.00  A,  Im2  =  2.00  A  case  performs  the  least  efficiently.  Factoring  both 
the  current  utilization  and  beam  divergence  efficiencies,  the  best  performing  configurations  were  the  Im \  = 
1.75  A,  Jm2  =  1.75  A  (least  divergence)  and  the  Im i  =  4.00  A,  Jm2  =  4.00  A  (best  current  utilization)  cases. 

The  krypton  power-matched  data  is  shown  in  Fig.  7.  For  krypton,  the  best  beam  focusing  was  produced 
by  the  4>div  optimized  magnet  current  setting  in  both  the  power-matched  and  volumetric  flow  rate  cases. 
The  inverse  relationship  between  divergence  efficiency  and  total  beam  current  observed  for  xenon  was  also 
seen  in  the  krypton  data.  While  the  Im i  =  4.00  A,  Jm2  =  2.00  A  case  produced  the  least  current  utilization 
efficiency  in  xenon,  it  provided  the  highest  value  in  the  case  of  krypton.  Trends  in  the  data  were  similar 
between  the  flow-matched  and  power-matched  cases.  In  general,  the  power-matched  cases  produced  higher 
efficiencies  compared  to  the  flow-matched  cases  suggesting  that  further  increases  in  mass  flow  rate  may 
improve  krypton  performance. 
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The  xenon  and  krypton  data  are  compared  in  Fig.  8.  Xenon  is  observed  to  have  about  a  17  %  divergence 
efficiency  advantage  compared  to  krypton  and  also  better  propellant  mass  utilization  efficiency.  The  higher 
beam  currents  for  xenon  also  suggest  that  xenon  has  a  better  propellant  utilization  efficiency  (which  factors 
into  thrust  and  specific  impulse  calculations).  However,  due  to  a  lack  of  charge  species  population  data,  the 
propellant  mass  utilization  advantage  cannot  be  quantified  precisely. 


Table  2.  Faraday  probe  derived  quantities. 

Im  Config. 


I ml 

(A) 

I m2 
(A) 

Prop. 

Match 

rna 

(mg/s) 

P 

(W) 

la 

(A) 

h 

(A) 

Jb 

Ia 

Jb 

la 

Rel.  to  Xe 

<t>DIV 

< pDIV 

Rel.  to  Xe 

<0  >j 

deg. 

1.07 

2.45 

Xe 

2.50 

582 

1.94 

1.833 

0.945 

0.795 

26.9 

0.56 

2.16 

Kr 

Va 

1.59 

495 

1.65 

1.419 

0.860 

-9% 

0.653 

-18% 

36.1 

0.56 

2.16 

Kr 

P 

1.81 

582 

1.94 

1.706 

0.879 

-7% 

0.671 

-16% 

35.0 

Xe 

2.50 

579 

1.93 

1.844 

0.955 

0.794 

27.0 

1.75 

1.75 

Kr 

Va 

1.59 

504 

1.68 

1.507 

0.897 

-6% 

0.645 

-19% 

36.6 

Kr 

P 

1.77 

579 

1.93 

1.738 

0.900 

-6% 

0.653 

-18% 

36.1 

Xe 

2.50 

609 

2.03 

1.877 

0.925 

0.765 

29.0 

4.00 

2.00 

Kr 

Va 

1.59 

504 

1.68 

1.573 

0.937 

+1% 

0.634 

-17% 

37.2 

Kr 

P 

1.84 

609 

2.03 

1.912 

0.942 

+2% 

0.641 

-16% 

36.8 

Xe 

2.50 

579 

1.93 

1.901 

0.985 

0.773 

28.4 

4.00 

4.00 

Kr 

Va 

1.59 

522 

1.74 

1.544 

0.887 

-10% 

0.649 

-16% 

36.3 

Kr 

P 

1.77 

579 

1.93 

1.794 

0.929 

-6% 

0.658 

-15% 

35.8 

Retarding  Potential  Analyzer 

Ion  energy-per-charge  (eV /q)  distribution  functions  (EDF’s)  were  measured  for  each  case  with  a  retarding 
potential  analyzer  at  radii  of  40  cm  and  80  cm.  Data  was  taken  in  5°  increments  from  9  =  -120°  to  9  =  120°. 
Voltage  on  the  ion  retarding  grid  was  swept  between  0  and  500  V  in  2  V  increments  relative  to  chamber 
ground.  The  difference  between  the  local  plasma  potential  of  the  measurement  location  and  chamber  was 
expected  to  be  less  than  1%  of  the  anode  voltage.  This  offset  was  neglected  in  energy  distribution  calculations. 
Each  current-voltage  trace  was  individually  fitted  with  a  smoothing  spline  before  being  differentiated  to 
calculate  the  ion  energy  distribution  function.  Each  EDF  was  area  normalized  to  unity  and  a  mean  energy 
was  calculated.  The  EDF  data  at  40  cm  matched  closely  with  80  cm  for  all  operating  conditions. 

Figure  9(a)  shows  the  EDF’s  for  xenon  propellant  at  a  radius  of  80  cm  for  various  angular  locations.  At  9 
=  0°,  a  primary  beam  ion  energy  peak  is  observed  near  the  anode  potential.  The  highest  energy  peak  is  for 
the  Im i  =  4.00  A,  /m2  =  4.00  A  case  at  290  eV  and  the  lowest  peak  occurs  at  268  eV  for  the  (j> div  optimized 
case.  As  theta  increases,  low  energy  features  appear  in  the  EDF’s.  At  9  =  45°  a  significantly  higher  fraction 
of  low  energy  ions  exist  for  the  Im \  =  1.75  A,  Im 2  =  1.75  A  case  and  the  4>div  optimized  case  than  for  the 
Im i  =  4.00  A,  Im2  =  2.00  A  case.  This  trend  seems  to  correlate  with  the  divergence  efficiency  data  from  the 
Faraday  probe  as  the  more  focused  configurations  produce  less  high  energy  ions  at  high  angles.  At  9  =  75° 
and  9  =  90°  the  EDF’s  are  dominated  by  low  energy  peaks.  However,  a  small  fraction  of  primary  ions  is 
still  observed  near  anode  potential.  Figure  9(b)  displays  the  trend  of  the  EDF  mean  over  angular  position. 
At  high  angles,  the  relationship  between  plume  divergence  and  distribution  mean  is  apparent.  The  more 
divergent  beams  exhibit  higher  mean  energy. 

Krypton  energy  distribution  functions  are  shown  in  Fig.  10(a).  For  krypton,  the  spread  of  energy  peaks 
among  the  different  cases  at  9  =  0°  was  smaller  than  for  xenon.  Opposite  of  xenon,  the  4>div  optimized 
case  had  the  highest  energy  peak  at  280  eV.  The  lowest  peak  was  at  270  eV,  occurring  for  the  Im\  =  4.00  A, 
Im2  =  4.00  A  case.  Mean  energy  data  in  Fig.  10(b)  shows  the  same  relationship  between  plume  divergence 
and  mean  energy  distribution  observed  for  xenon  also  applies  to  krypton. 

Krypton  and  xenon  EDF’s  are  plotted  together  in  Fig.  11.  One  of  the  most  noticeable  differences  between 
xenon  and  krypton  EDF’s  is  that  krypton  has  a  narrower  primary  energy  peak  indicating  a  smaller  spread 
in  primary  ion  energy.  Also  due  to  the  higher  divergence  of  krypton,  a  larger  relative  population  of  high 
energy  ions  is  found  at  high  angles  as  seen  in  the  9  =  60°  and  9  =  75°  EDF  plots. 

Voltage  utilization  efficiency  was  calculated  using  Eq.  5,  where  the  EDF  distribution  mean  was  ion 
acceleration  voltage.  The  voltage  utilization  efficiency  was  calculated  twice  for  each  operating  condition 
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using  Faraday  probe  and  RPA  data  from  both  40  cm  and  80  cm.  These  two  separately  calculated  values 
provide  an  indication  of  repeatability  of  this  technique  for  calculating  voltage  utilization  efficiency  from 
probe  data. 

The  approach  taken  to  evaluate  voltage  utilization  efficiency  in  this  study  is  different  than  the  more 
common  method  of  only  using  the  ratio  of  most  probable  energy  to  anode  voltage.3, 9- 10  Other  researchers 
use  the  most  probable  energy  approach  because  probe  noise  uncertainty  results  in  inaccurate  integration  of 
the  energy  distribution  function.9  Therefore,  the  calculated  mean  energy  is  less  representative  of  the  voltage 
utilization  than  the  most  probable  energy.3  However,  RPA  measurements  in  this  study  were  found  to  be 
highly  repeatable,  similar  for  both  40  cm  and  80  cm  radius  at  a  constant  angular  location,  and  symmetrical 
about  0  =  0°.  Although,  systematic  uncertainty  exists  in  these  measurements,  uncertainty  from  random 
noise  appeared  to  be  low.  In  this  study,  voltage  utilization  is  intended  to  be  used  as  a  comparison  metric 
among  cases  rather  than  a  figure  used  in  calculating  total  anode  efficiency. 

Table  3  contains  the  calculated  values  of  voltage  utilization  efficiency  for  each  case.  Calculations  for  the 
different  radii  were  highly  repeatable  with  less  than  2%  deviation.  For  xenon,  the  best  voltage  utilization 
efficiency  was  calculated  for  the  Iml  =  4.00  A,  Im2  =  4.00  A  case.  In  general,  power-matched  krypton  cases 
produced  higher  voltage  utilization  efficiency  compared  to  the  flow-matched  cases.  The  Im i  =  4.00  A,  Im2 
=  2.00  A  case  produced  the  best  voltage  utilization  efficiency  for  krypton.  Voltage  utilization  efficiency 
apperared  to  be  higher  for  cases  with  higher  beam  current. 


Table  3.  RPA  Derived  Quantities 

Im  Config. 


1ml 

(A) 

I  m2 

(A) 

Prop. 

Match 

rna 

(mg/s) 

P 

(W) 

la 

(A) 

AV 

Va 

(40  cm) 

rel. 

to  Xe 

AV 

Va. 

(80  cm) 

rel. 

to  Xe 

1.07 

2.45 

Xe 

2.50 

582 

1.94 

0.761 

0.775 

0.56 

2.16 

Kr 

Va 

1.59 

495 

1.65 

0.785 

3% 

0.782 

1% 

0.56 

2.16 

Kr 

P 

1.81 

582 

1.94 

0.799 

5% 

0.779 

1% 

Xe 

2.50 

579 

1.93 

0.814 

0.819 

1.75 

1.75 

Kr 

va 

1.59 

504 

1.68 

0.743 

-9% 

0.749 

-8% 

Kr 

P 

1.77 

579 

1.93 

0.806 

-1% 

0.806 

-2% 

Xe 

2.50 

609 

2.03 

0.784 

0.793 

4.00 

2.00 

Kr 

Va 

1.59 

504 

1.68 

0.741 

-5% 

0.744 

-6% 

Kr 

P 

1.84 

609 

2.03 

0.814 

4% 

0.806 

2% 

Xe 

2.50 

579 

1.93 

0.836 

0.823 

4.00 

4.00 

Kr 

Va 

1.59 

522 

1.74 

0.737 

-12% 

0.737 

-10% 

Kr 

P 

1.77 

579 

1.93 

0.799 

-4% 

0.783 

-5% 

Discharge  Current  Oscillation  Effects 

Anode  current  data  was  measured  in  the  time  domain  using  an  oscilloscope  and  in  the  frequency  domain 
using  a  signal  analyzer  for  all  tested  operating  conditions.  Magnet  current  settings  were  observed  to  strongly 
influence  oscillatory  behavior  of  anode  current.  Figure  12  shows  the  anode  current  behavior  for  xenon.  Both 
the  7mi  =  1.75  A,  Im2  =  1.75  A  and  the  <!>div  optimized  case  were  observed  to  have  a  highly  oscillatory 
nature  and  harmonic  frequencies.  The  most  quiescent  case  was  the  Im i  =  4.00  A,  Im2  =  4.00  A.  The  peak 
frequency  was  highest  for  the  quiescent  case  (41.1  kHz)  and  decreased  as  the  oscillation  magnitude  increased 
(37.2  kHz  for  the  4>div  optimized  case). 

Figure  13  shows  the  anode  current  oscillations  for  the  krypton  power-matched  cases  and  Fig.  14  displays 
this  same  data  along  with  the  xenon  data  for  comparison.  Overall,  the  krypton  oscillation  magnitudes  were 
much  smaller  than  for  xenon.  Similar  to  xenon,  the  krypton  4>div  optimized  case  produced  the  largest 
magnitude  oscillations.  The  Im\  =  4.00  A,  Im2  =  2.00  A  case  was  most  quiescent.  However,  in  contrast 
to  xenon  the  most  quiescent  krypton  cases  had  the  lowest  peak  frequency  (20.5  kHz)  and  peak  frequency 
increased  with  increasing  oscillation  magnitude  (41.6  kHz  for  the  divergence  optimized  case). 
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Table  4  summarizes  the  anode  current  data  for  all  the  operating  conditions.  The  beam  current  appears 
to  be  strongly  related  to  the  anode  current  oscillations.  For  xenon,  high  beam  current  cases  tend  to  have 
low  oscillation  magnitudes  and  high  peak  frequency.  Whereas  for  krypton,  high  beam  current  cases  tend  to 
have  both  low  oscillation  magnitude  and  low  peak  frequency.  Consequently,  low  divergence  cases  (having 
low  beam  current)  produce  large  anode  oscillations. 


Table  4.  Anode  Current  Data 

Im  Config. 


I ml 

(A) 

I m2 
(A) 

Prop. 

Match 

rha 

(mg/s) 

P 

(W) 

la 

(A) 

h 

(A) 

Jb 

Ia 

<t>DIV 

AV 

Va 

(80  cm) 

<?Ia 

(A) 

/peak 

(kHz) 

Pi* 

(dB-A) 

1.07 

2.45 

Xe 

2.50 

582 

1.94 

1.833 

0.945 

0.795 

0.775 

0.78 

37.2 

-18.9 

0.56 

2.16 

Kr 

Va 

1.59 

495 

1.65 

1.419 

0.860 

0.653 

0.782 

0.23 

37.5 

-34.6 

0.56 

2.16 

Kr 

P 

1.81 

582 

1.94 

1.706 

0.879 

0.671 

0.779 

0.18 

41.6 

-43.4 

Xe 

2.50 

579 

1.93 

1.844 

0.955 

0.794 

0.819 

0.80 

36.4 

-18.1 

1.75 

1.75 

Kr 

Va 
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Conclusion 

This  investigation  compared  plume  data  from  a  medium  power  Hall  effect  thruster  operating  with  xenon 
and  krypton  propellant.  Faraday  probe  data  was  used  to  calculate  divergence  efficiency,  beam  current,  and 
beam  current  utilization  efficiency.  Energy-per-charge  distribution  functions,  measured  with  an  RPA,  were 
used  to  calculate  voltage  utilization  efficiency.  Anode  current  data  were  analyzed  in  the  frequency  domain 
to  correlate  oscillatory  behavior  with  beam  properties. 

The  methodology  of  this  plume  study  was  different  than  typical  performance  analysis  studies  in  two 
ways.  First,  probe  data  was  taken  at  multiple  radii.  Faraday  probe  data  were  typically  taken  between  7  (for 
Kr  testing)  and  13  (for  Xe  testing)  different  radii.  Divergence  efficiency  was  observed  to  exhibit  a  strong 
linear  relationship  with  probe  sweep  radius.  Retarding  potential  analyzer  data  was  taken  at  two  radii  for 
each  operating  condition.  Energy-per-charge  distribution  functions  remained  similar  at  both  radii. 

The  second  difference  was  that  charge-weighted  integration  of  mean  energy  was  used  to  calculate  voltage 
utilization  efficiency  instead  of  using  only  the  most  probable  energy  from  EDF’s  measured  near  0  =  0°. 
This  method  of  calculation  is  more  analytically  rigorous,  but  is  possibly  subject  to  greater  uncertainty  from 
integrating  EDF’s  that  have  higher  uncertainties  at  low  energy.  A  more  detailed  study  on  the  accuracy  of 
RPA  derived  EDF’s  for  a  medium  power  Hall  thruster  will  have  to  be  conducted  to  determine  which  voltage 
utilization  efficiency  calculation  technique  is  most  appropriate. 

Faraday  probe  data  revealed  that  krypton  has  a  divergence  efficiency  value  approximately  15-20%  less 
than  xenon.  It  also  showed  that  krypton  current  utilization  efficiency  is  5-10%  less  than  that  of  xenon. 
Different  radial  magnetic  field  strength  configurations  had  a  significant  effect  on  diverence  and  current 
efficiency  values. 

Differences  between  voltage  utilization  efficiency  for  xenon  and  krypton  were  small  (approximately  2% 
less  for  krypton).  Voltage  utilization  efficiency  was  observed  to  be  linked  to  beam  current,  where  cases  with 
higher  beam  current  had  higher  voltage  utilization  efficiency.  Retarding  potential  analyzer  data  showed 
that  significant  populations  of  high  energy  ions  existed  at  high  angles  for  operating  conditions  with  high 
divergence. 

Different  magnetic  fields  produced  variations  in  anode  current  oscillatory  behavior.  Anode  current  oscil¬ 
lation  magnitude  was  observed  to  be  significantly  smaller  for  krypton  than  xenon.  Operating  conditions  with 
high  beam  current  tended  to  produce  larger  magnitude  oscillations.  For  xenon,  peak  frequency  decreased 
with  oscillation  magnitude,  but  for  krypton  the  opposite  effect  was  observed. 
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Besides  the  efficiency  losses  studied  here  (relating  to  anode  efficiency),  the  propellant  mass  utilization 
efficiency  (rhi/rha)  appeared  to  be  significantly  less  for  krypton.  The  propellant  mass  utilization  efficiency 
cannot  be  calculated  accurately  without  charge  species  fractions  and  was  not  quantified  due  to  the  lack  of 
E  x  B  data  in  this  study.  Specific  impulse  is  directly  proportional  to  mass  utilization  efficiency.  Losses  in 
mass  utilization  efficiency  would  directly  counteract  the  prized  high  ion  velocity  of  krypton  in  the  calculation 
of  specific  impulse.  Considering  this  loss  combined  with  the  higher  beam  divergence  loss,  the  specific  impulse 
for  krypton  may  not  be  better  for  krypton  than  xenon  for  this  medium  power  Hall  thruster. 

Ultimately,  a  more  thorough  experimental  study  that  includes  charge  species  measurements  from  an 
E  x  B  probe  and  thrust  data  from  a  thrust  stand  is  needed  to  make  more  complete  and  accurate  evaluations 
of  the  anode  efficiency  and  thruster  performance.  Also,  a  more  rigorous  survey  of  operating  conditions  that 
includes  increasing  anode  flow  rate  and  discharge  potential  should  be  performed  to  find  operating  conditions 
that  maximize  the  efficiency  and  performance  of  krypton.  These  conditions  would  provide  more  relevant 
comparisons  to  the  established  xenon  operating  conditions.  Initial  observations  for  a  limited  number  of 
operating  cases  indicate  that  krypton  is  significantly  less  efficient  than  xenon  and  likely  does  not  offer  a 
significant  increase  in  specific  impulse  for  the  BHT-HD-600  Hall  thruster. 

Appendix:  Experimental  Data  Figures 
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(c)  Divergence  efficiency  calculated  at  each  probe  sweep  radius,  (d)  Beam  current  (integration  limit  90  deg.)  calculated  at  each 

probe  sweep  radius. 
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(e)  Interpolated  plume  charge  flux  contour. 


Figure  6.  Xenon  Magnet  Current  Comparison:  Faraday  probe  data  for  thruster  operating  at  \ '  —  .'1(H)  V,  and 
ma= 25.5  seem. 
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(c)  Divergence  efficiency  calculated  at  each  probe  sweep  radius,  (d)  Beam  current  (integration  limit  90  deg.)  calculated  at  each 

probe  sweep  radius. 


(e)  Interpolated  plume  charge  flux  contour. 


Figure  7.  Krypton  (power  matched)  Magnet  Current  Comparison:  Faraday  probe  data  for  thruster  operating  at 
Vd=300  V. 
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(a)  Charge  flux  vs.  theta,  r  =  60  cm. 


(c)  Divergence  efficiency  calculated  at  each  probe  sweep  radius. 


(d)  Beam  current  (integration  limit  90  deg.)  calculated  at  each 
probe  sweep  radius. 


(e)  Interpolated  plume  charge  flux  contour. 


Figure  8.  Xenon  vs.  Krypton  (power  matched)  Magnet  Current  Comparison:  Faraday  probe  data  for  thruster 
operating  at  V^=300  V. 
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-  Im1  =  1.75  A,  lm2  =  1.75  A 

lm1  =  1-07  A,  lm2  =  2.45  A 
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-  Imi  =  4.00  A,  lm2  =  4.00  A 
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(b)  Energy  Distribution  Function  Mean  Values,  r  =  80  cm 

Figure  9.  Xenon  Magnet  Current  Comparison :  RPA  data  for  thruster  operating  at  V^=300  V,  and  ma= 25.5  seem. 
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'mi  =0-56  A,  lm2  =  2.16  A 
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(b)  Energy  Distribution  Function  Mean  Values,  r  =  80  cm 

Figure  10.  Krypton  (power  matched)  Magnet  Current  Comparison:  RPA  data  for  thruster  operating  at  V^=300  V. 
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Figure  11.  Xenon  vs.  Krypton  (power  matched)  Magnet  Current  Comparison:  RPA  data  for  thruster  operating 
at  Vd=300  V. 
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Figure  12.  Xenon  Magnet  Current  Comparison :  Anode  current  data  for  thruster  operating  at  V^=300  V,  and  rha=25.5  seem. 
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Figure  13.  Kryptonfpower  matched)  Magnet  Current  Comparison:  Anode  current  data  for  thruster  operating  at  Ud=300  V. 
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Figure  14.  Xenon  vs.  Krypton  (power  matched)  Magnet  Current  Comparison:  Anode  current  data  for  thruster  operating  at  Ud=300  V. 
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